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The problem of cosmological production of gravitational waves is discussed in the framework of 
an expanding, spatially homogeneous and isotropic FRW type Universe with time-evolving vacuum 
energy density. The gravitational wave equation is established and its modified time-dependent part 
is analytically resolved for different epochs in the case of a flat geometry. Unlike the standard ACDM 
cosmology (no interacting vacuum), we show that gravitational waves are produced in the radiation 
era even in the context of general relativity. We also show that for all values of the free parameter, 
the high frequency modes are damped out even faster than in the standard cosmology both in 
the radiation and matter-vacuum dominated epoch. The formation of the stochastic background 
of gravitons and the remnant power spectrum generated at different cosmological eras are also 
explicitly evaluated. It is argued that measurements of the CMB polarization (B-modes) and its 
comparison with the rigid ACDM model plus the inflationary paradigm may become a crucial test 
for dynamical dark energy models in the near future. 


I. INTRODUCTION 

The recent direct detection of gravitational waves (GWs) by the LIGO collaboration [1] is a turning point in 
the field opening the new era of GW astronomy. This result provoked both surprise and excitement not only on 
the researchers working with the generation of GWs from astronomical objects but also of cosmological origin. 
The implications are promising because the weak interaction of the GWs with matter will provide information 
previously unattainable like the physics of strong relativistic gravitational fields (black holes coalescence), spe¬ 
cific signatures from the very early universe (much before CMB), and even tests of gravitational theories [2]. 
Potentially, cosmological GWs will furnish new and powerful key for accessing the early universe physics and 
their current problems like the mystery related with the nature of the unknown components in the Universe 
(dark matter and dark energy). 

In this concern, the simplest explanation for the present accelerating stage of the observed Universe is the 
existence of a dark energy (DE) component, in addition to the cold dark matter (CDM) [3-5]. Observationally, 
the most accepted candidate for DE is a rigid cosmological constant A with energy density p vac = A/87rG. 
At the level of general relativity (GR) and quantum held theory (QFT) in curved spacetimes, the A-terrn can 
be interpreted as a decoupled relativistic simple fluid (vacuum medium) obeying the equation of state (EoS), 
p v ac = -p v ac = constant [6]. 

The recent Planck2015 results confirmed that the cosmic concordance model (ACDM) - a hat Universe filled 
with cold dark matter (CDM), baryons and a A-terrn - is in good agreement with the currently available 
cosmological observations [7]. However, from the theoretical viewpoint, the unsettled situation in the particle 
physics/cosmology interface in which the cosmological upper bound (p vac < lCU 47 GeU 4 ) differs from naive 
theoretical QFT expectations (p vac ~ 10 71 GeV 4 ) by more than 100 orders of magnitude, originates an extreme 
hne-tuning problem, the so called cosmological constant problem [6, 8]. Another unsolved mystery from hrst 
principles is why the vacuum energy density is so close to the matter energy density which is sometimes referred 
to as cosmic coincidence problem [9-11]. 

In order to alleviate the cosmological constant and coincidence problems, and also some mild tensions of the 
ACDM model [12-15] many authors have proposed dynamical DE models (for a recent discussion see [16]). In 
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general grounds, these models can also be implemented trough a time-evolving vacuum energy density coupled 
with the remaining cosmic components. Theoretically, such models are suggested by the general form of the 
effective action in QFT in curved spacetimes [17] . However, there are also some attempts to represent such A 
models by a scalar Held [18, 19], a Lagrangian description based on the so-called F(R,T ) gravity where R is 
the curvature scalar and T is the trace of the energy-momentum tensor [20, 21]. In such frameworks, different 
phenomenological decay laws for a time-dependent A have also been proposed and their predictions confronted 
with the available observational data [22-33]. 

Broadly speaking, the basic idea of a time-evolving A-term is quite simple. The vacuum EoS implies that 
Pvac is constant only when the vacuum component is separately conserved. This is not the case when an 
interacting mixture is considered as in the present work. To be more specific, consider the corresponding 
Einstein field equations for a system formed by the interacting mixture (matter plus the vacuum fluid) with 
energy momentum tensors and respectively. Since the Einstein tensor is divergenceless, the 

total energy conservation law in the direction of the four-velocity, u^(T^f at ) + T (vac));- = °> reads: 

M M T (mat);^ = P + 3H(p + p) = -p vac = - ^ , (1) 

where a FRW geometry was assumed and in the last equality the vacuum EoS was used (see (5)). Note that 
the energy of the material component is not conserved independently. The running vacuum component is 
continuously transferring energy to the fluid component. Of course, when the energy density of the vacuum 
medium becomes constant it does not contribute to the above balance equation for the material component. 
In other words, if A (t) is the constant (bare) A(,, the energy of the material medium is separately conserved 
thereby recovering the standard result. 

On the other hand, the recent LIGO results are pressing the search on the production of relic GWs (tensor 
modes) both at very low and high frequencies. In the observational front, the major challenge is the detection 
of the B-modes polarization anisotropies in the cosmic background radiation (CMB) provoked by tensor fluctu¬ 
ations produced in the inflationary regime. The expected level of such anisotropies is small but presumably it 
will be accessed by the ongoing (BICEP and Keck Array experiments, South Pole Telescope, etc.) and future 
probes through a new generation of high sensitivity instruments. One of them is the QU bolometric interferom¬ 
eter for cosmology (QUIBIC) whose basic goal is to measure the B-modes at angular intermediate scales based 
on the innovative technological concept of bolometric interferometry (see, for instance, [34]). In principle, the 
presence of a time-evolving A-term must effect the B-modes in comparison with the traditional ACDM plus the 
inflationary paradigm. As a consequence, the influence of a dynamical vacuum on the adiabatic amplification 
of GWs is a possibility that needs to be investigated. This is one of the main aims of the present paper. 

In this context we investigate here the production of primordial GWs in the framework of a simple A (H) 
decaying vacuum cosmology that recovers the ACDM model as a particular case. The energy density and power 
spectrum are analytically derived. Interestingly, unlike the cosmic concordance model (ACDM), there exists 
adiabatic amplification of GWs during the radiation epoch. As we shall see, such a result holds for generic 
decaying vacuum models and must remain valid for all analyzes based on realistic running vacuum cosmologies. 


II. THE MODEL: BASIC EQUATIONS 

Let us now consider that the Universe is well described by a flat FRW geometry. In the co-moving coordinate 
system, the background fine element reads (c = 1): 


ds 2 = - dt 2 + a 2 (t) dl 2 = a 2 (if) (-dr? 2 + dl 2 ) , (2) 

where the 3-space metric is dl 2 = dx 2 + dy 2 + dz 2 = dijdx l dx\ i,j,k = 1,2,3 and t, rj are, respectively, the 
cosmic and conformal times, related by 


dt = a drj . 


( 3 ) 


The Einstein field equations for a two-fluid mixture (matter + vacuum) is given by: 
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W v -\r g^ = 87 tG( TZt + TZ) , (4) 

with and R , being the Ricci tensor and Ricci scalar, respectively. The energy-momentum tensors (matter 
and vacuum) read: 


TZt = (p + p) / - p , TZ = p vac sT , (5) 

where p , p, u ^ are, respectively, the energy density, pressure and four-velocity of the fluid, whereas p vac = 
A(t)/87rG, is the vacuum energy density. The field equations in the above background take the form [23, 25, 2' 

a 2 

8nG p + A (t) = 3— , (6) 

a z 

• 2 

Sir Gp - A (t) = -2-—°— (7) 

a or 

where a dot means derivative with respect to the cosmic time t. In order to solve the above equation we need 
to know both the functional form of A(f) and the EoS. 

Many phenomenological functional forms have been proposed in the literature for describing a time-varying 
A(<) vacuum [27]. Based on dimensional arguments, Carvalho et al. [23] shown that a natural dependence is 
A oc H 2 . Later on, Waga added the bare cosmological constant to this H 2 dependence studying some physical 
consequences of law A {H) = At, + (3H 2 , where (3 is a dimensionless free parameter [24]. Finally, such a functional 
dependence was derived within a renormalization group approach based on QFT in curved spacetime by Sola 
and Shapiro [17]. 

In what follows, in order to discuss analytically some subtleties related to the decaying vacuum contribution 
for the stochastic background of GWs, it will be assumed here that the phenomenological A(£7)-term is the one 
adopted by the above authors [17, 23, 24] which is also a particular case of the class studied by the authors of 
Refs. [29, 30]: 


A (H) = A b + 3 pH 2 , 


where H = a/a is the Hubble parameter, and the factor 3 was added for mathematical convenience. 
By assuming the w-law EoS 


( 8 ) 


P = up, (9) 

where w is a different constant for each era, and combining Eqs. (6)-(9), we find that the scale factor is driven 
by the following differential equation [35]: 


ad + A a 2 — ^ ^ ^ A b a 2 = 0 , (10) 

where we have introduced the convenient short notation 

a = 3(1-M(1 — a — 2 (11) 

The case j3 = 0 reduces to the standard ACDM equation for different eras. In the absence of A b , the free 
/3 parameter was first constrained long back ago based on big-bang nucleosynthesis studies to be j3 < 0.16 
36, 37]. Actually, it is even smaller when the constant vacuum term is also considered in the same decaying 
law (/3 ~ 10" 3 ) [32], 

In terms of the conformal time ( 77 ), the scale factor equation (10) can be rewritten as (see [38] for /3 = 0): 
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a" a' 2 (1 + w) 2 

— + (A - 1 )— - - A b a 2 = 0 . 

a a z 2 


( 12 ) 


where primes denote derivatives with respect to 77 . The last term in the above equation is negligible at early 
times. Actually this happens for redshifts of the order of few. In this case, the general solution for the scale 
factor can be written as: 


a{rj) = Ci(A ?7 - c 2 ) 1/A , (13) 

where c\ and C 2 are integration constants. As usual, the solution for each era is specified by the corresponding 
value of w, namely: inflation (wi n f = —1), radiation-vacuum (w ra d = 1/3) and matter-vacuum (w mat = 0). For 
each era we have also assumed that the vacuum decays only on the dominant component. Note also that the 
early exponential inflation is not modified by the decaying vacuum since w = — 1 implies A = — 1 so that d > 0 
regardless of the values of f3 (see Eqs. (10) and (11)). 

In order to find the integration constants for each era we must use the continuity junction conditions for the 
transition times between each era, a n (r)i) = a n +i(jii) and a' n {rji) = a' n+1 (r]i). With this we have: 


—hv 1 V < 0, V < m 


a{rf) = < 


^^Or(A ra d77 77rad) ^ rad 5 


Vi < V < Veq 


ka 0m (A ma , t r) - Tjryi) 1 / Amat , 77 > 7? eq 


(14) 


where U is a constant and the parameter A Q corresponds to the value oj a . The transition time between inflation 
and radiation era is 771 , and between radiation and matter is 77 eq . The values of the integration constants are 


Vr&d (A ra d T l)77l, 
aor = (-r/i) _(1+1/Arad) , 

^?mat = (^mat ^rad)^7eq “1“ Vr&di 

(^rad^?eq ^7rad) ^ rad 

^Om — ^rad ~T~\ \i / A 

(^mat^eq ^/mat) ' mat 


(15) 


Now we have a complete solution for the scale factor which depends on the parameter f3. In the special case 
of no decaying vacuum (/3 = 0), the quantity A as given by (11) reduces to the standard definition and all the 
expressions appearing in Ref. [39] are recovered. For future calculations it is important to estimate values of 
the transition times 771 and 77 eq . In order to do that, let us compare the scale factors of different cosmological 
eras. Comparing the scale factor at the present time 770 with the scale factor at the end of inflation we have 
a(r]o)/a(r/i) ~ 10 21 and for the end of the radiation era 0 ( 770 )/ a(r} eq ) ~ 10 4 [39]. Using the solution for the 
scale factor (14) and solving the equations system we obtain that 77 ! ~ —10 -17 and ? 7 eq ~3x 10 -3 . Although 
representing a crude approximation since the value of rj eq must depend on the /3 parameter, throughout this 
paper we adopt this values. 

In the Figure 1 we show schematically the behavior of the scale factor for some selected values of (3. As 
expected, for a given value of ui, the expansion grows faster for higher values of /?. This happens because the 
vacuum component contributes with a negative pressure. 


III. COSMOLOGICAL TENSOR PERTURBATIONS 

In the conformal time, a classical tensor metric perturbation in the FRW flat geometry given by Eq.(2) can 
be written as: 


ds 2 = a 2 (ri)[—dr] 2 + (t>^ + h i j)dx l dx :i } 


(16) 
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FIG. 1: Evolution of the scale factor for different regimes and some selected values of /3 as indicated in the figure. The 
expansion rate is faster for higher values of /3. 


where the perturbation h l] is small, h l3 1 <C 1, and satisfy the well known [10] symmetry (transverse-traceless) 
and gauge constraints, namely: h 03l = 0, h\ = 0, V J h l3 = 0. 

In the context of GR, the evolution of the tensor perturbations is given by the standard form: 


+ - V 2 ^' =0. (17) 

This happens because both components in the mixture (matter plus vacuum) have the perfect fluid isotropic 
form [41, 42]. The general solution can be expressed in terms of Fourier expansion: 


hijiv,*-) = 


\Zl6nG 
(2tt) 3 / 2 , 




E ; 

r=+, x 


;( n ) 


h n (r])e in ' x c n + h* n {i))e~ 


(18) 


where h n {rf) are the mode functions, n is the comoving wave vector, c n and c„ are complex numbers. 
The polarization tensor, n), is symmetric (eij(n) = ejj(n)), traceless (eii(n) = 0), and transverse 
( niEij(n ) = 0). We also choose a circular-polarization basis in which e l3 (n) = n ))*, and normalize the 

basis Y,i,j Gi( n )(Gj( n ))* = 2 ^ s - 

The comoving wave number n = |n| is related with the physical wave number k by 


n= |n| = =ka(r,). (19) 

A 

Now, by inserting the solution (18) in (17) it is readily seen that the temporal part decouples thereby giving 
the evolution equation for the conformal time modes: 


h n (v)" + 2 — hnivY +n 2 h n (r]) = 0. (20) 

a 

r r 

Using the auxiliary function = h n (r])a(r]) the above equation assumes the first obtained by Grishchuck 

[39] which is independently satisfied for each polarization r = +, x: 


r // 



( 21 ) 


Therefore, by assuming that the vacuum component is smooth we see that the standard wave equation is not 
modified. Given the solutions for the scale factor a(r/) in different eras, we can solve (21) for each mode n. 
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H >?cq 

V 

FIG. 2: Potential, V(rj) = a”/a for some values of /3. Note that there is no adiabatic amplification for /3 = 0 and uj = 1/3 
(standard radiation era) since in this case a"/a = 0 (see also Eqs. ( 11 )-( 12) ). 


It represents an harmonic oscillator with variable frequency determined by the evolution of the Universe and 
describes different behaviors for the high and low frequency regimes (with and without a vacuum component). 
Once the solutions for p(if) for the different cosmic eras has been calculated, it is immediate to obtain the 
associated physical quantities like the wave amplitude, energy density and power spectrum. 

An important quantity driving the behavior of the primordial GWs is the “potential” V(rj ) = a"/a appearing 
in equation ( 21 ) (it should be recalled that the name “potential” come from the mathematical analogue with 
the stationary Schrodinger equation). The relation between the potential and the wave-number determines the 
behavior of the limit solutions for p(rj). 

For the times when n 2 \V\ holds, the solution of (21) is oscillatory, p oc e ±m7) , so that the high-frequency 
waves are diluted by the cosmic expansion h = e ±lnr, /a. In the opposite limit, n 2 <C \V\, we have p ex a, and, 
consequently, the low-frequency waves obey h = constant. The effect of the potential is to avoid the damping 
of the waves due the universe expansion. The net effect is that the perturbations are relatively enhanced, a 
phenomenon commonly referred to as adiabatic amplification [43]. Note that in the limit case (/3 = 0, ui = 1/3), 
that is, in the standard radiation phase, we see that A = 1 and, therefore a" = 0. It thus follows that in the 
radiation era (771 < p < rj eq ), the potential vanishes identically (V(p) = 0). Physically, this means that there is 
no adiabatic amplification of GWs during the standard radiation phase. 

In Figure 2 we show the behavior of the potential for the different eras. As it will be discussed next, for 
ui = 1/3 and /? 7 ^ 0 GWs are produced so that low frequency modes can be slightly amplified even during the 
radiation phase. Let us now discuss the solutions of the wave equation for the different eras. 


IV. GRAVITATIONAL WAVE SOLUTIONS FOR DIFFERENT ERAS 

A. Inflation Era 

In the particular case of cu = cu; n f = —1 we have an exponential inflation with a potential, a"/a = —2/p 2 , 
regardless of the values assumed by the /? parameter. This case has already been studied in [44]. The scale 
factor a,i n f( 77 ) = —li/r) gives a positive constant Hubble parameter Hj = l~ . In the inflation era the GW 
equation ( 21 ) can be written as: 


Minf "F ( n ^ ] Pini — 0, 


( 22 ) 


where for simplicity, we have suppressed the polarization index r. The general solution of the above equation 
can be expressed in terms of Bessel’s functions 


Pini(n,r]) = y/r\ [A, J_ 3 / 2 (w?) + BiJ 3/2 (np)}. 


(23) 
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We have to specify some conditions to calculate the integration constants Aj and B j. In the inflation era, 
the limit for high frequencies must reach the so-called adiabatic vacuum linin^oo /z oc e~ lnr> [45]. Using this 
condition and also in its first derivative the constants reduce to A,; = isJir/2 and Bi = — yjir/2. After doing 
some algebra we have the normalized (/q n f /z[* f — /n* n{ lAnf = *) solution 


rj) 



(24) 


Knowing the full expressions for /i and a is easy to calculate the power V and energy spectrum fl gw . For details 
see the Appendix B. 


V in i(n, rj) 

fi&V*7) 


16G 2 2 \ 

rif * 1 + B ” 


8 G 4 4 


3nh 


n V 2 + 


n 2 rj 2 


(25) 


These are standard results already studied, for details see [42]. In particular for long wavelengths X^> k = Hj 1 
the power spectrum is flat and proportional to Hi , V\ n i = 160 H 2 . For our purposes the calculation of /q n f is 
important to obtain a complete solution for /r rad in the radiation era that will be shown in the next section. 


B. Radiation Era 

In the radiation era we have u = w rad = 1/3 and A ra( j = 1 — 2/3. From the second solution given by Eq. (14), 
it is readily seen that the Eq. (21) now takes the form: 


/Gad A 



1 - A 


rad 


[A rad (r/ - J 7 i) - rj\Y 


Mrad — 0 • 


(26) 


At this point, we have to stress the first important result of the present paper: The potential for this case is 

V = —2(3(a'/a) 2 ^ 0. In the particular case of /? = 0 (no decaying vacuum) we obtain the well-known result 

V = 0 of no GW amplification in the radiation era. In the decaying vacuum models this does not hold anymore, 
the V/ 0 condition implies that always the radiation will contribute to the primordial GW spectrum today. 

The general solution of the last equation is: 


/irad (n,T/) = yj A rad 7y - 7/i(A rad + 1) 


(27) 


jy J f~y 


*rad 


(A rad ?/ - r/i(A rad + 1)) + B r J c 


n 


A 


rad 


(^rad?7 ^?l(^rad H” 1)) 


where a r = The continuity junction conditions, /q n f(? 7 i) = Hr&d(vi) and = At( ad (/ 7 i), must be 

used to calculate the integration constants A r and B r . Making the calculations we have: 


A r 

B r 


e~ m?? 1 7 rsec( 7 r/A rad ) 
A rad (—2n?7i ) 3 / 2 

e~ m?? 1 7 rsec( 7 r/A 1 . ad ) 

A rad (—2n?7i ) 3 / 2 


[nrji(i - to/i) Ji_ Qr . (77*) - ( 2 i - 77771 (2 + ini] 1) - A rad (i 
[-in 2 rjlJ aT (77*) (inr) - n 2 r]j)J ar+ i (77*)] . 


kr]i))J-oc r (V*)], 


where 77 * = — d/ 23 -. Note that the above solution for the radiation era is also normalized since it satisfies 
A*rad /4 ad — /x/ ad /7r ad = *• These expressions allow us to obtain the amplitude of the perturbations. 

Figures 3a and 3b display the evolution of the amplitude |/i rad | as a function of the conformal time 77 and 
some selected values of /3. As discussed above, the behavior of low frequency (Figure 3a) and high frequency 
(Figure 3b) modes are quite different since the former are amplified (even during the radiation phase) while 














FIG. 3: The amplitude of the GW in the radiation era. a) Modulus of the mode function |/i ra d| as a function of the 
conformal time for some selected values of fi and a fixed comoving low frequency, n = 10. b) The same plot of figure a) 
but now for a fixed high frequency, n = 10 5 . Note that in the low frequency regime decaying vacuum models amplify the 
perturbations since the amplitude is higher in comparison to the case /3 = 0. However, high frequency modes are always 
damped out regardless of the value of /3 (see main text), c) Modulus of the mode function |/i ra d| as a function of n and 
the same selected values of /3. 


the later are damped. Surprisingly, we see that the high frequencies modes are damped even faster than in the 
standard case (/3 = 0). 

How can such a result be understood? The basic point here is that in the high frequency limit the term 
a"/a can be neglected and h(rj) = a -1 (?y)/Lt(? 7 ). Therefore, since the solution of /r is an oscillating function, the 
amplitudes are damped out even more intensively (in comparison to /3 = 0) since the scale factor expands faster 
for higher values of (3 ^ 0 (see Figure 1). 

Conversely, in the low frequency regime we find exactly the opposite behavior. Indeed, due to the condition 
n 2 -C a"/a, the solution for low frequencies is /i oc a so that the perturbations, remains nearly constant (see 
Figure 3a). Is exactly in this regime (low frequencies) that the amplification occurs even during the radiation 
phase. Finally, in figure 3c we show the behavior of |/i ra d| as a function of the frequency for a fixed time. 


C. Matter Dominated era 

In the matter dominated phase, the EoS is to = w ma t = 0 which implies that A mat = (1 — 3/3)/2. In this case, 
the scale factor reduces to: 


^m a t(^?) — li 


&0 n 


1 — 3 p 


V ^7mat 


\ 2 /( 1 — 3 / 3 ) 


(28) 


1_« 2(l-ff) 1-ft 

where r? mat =- ^-i ) eq + 2(1 - and a 0m = (—i?i) ^ [(1 - 2/3)r) eq - 2(1 - 0)^] i-^+w* 2 . Calculating 

the potential and substituting in (21) we have 


Mma.t + 



2 + 6/3 


[(3/3 - 1)77 + 2?7 ma t] : 


Mm a t — 0? 


solving it we obtain the general solution 


(29) 


M mat 


— \/ZY 


mat V ^7mat 


A T 

rt 


i(A n 


A 


,tV ^7mat) ^ 


H - BmJ—o 


^(A ma t TJ ?7mat) 

Amat, 


(30) 


with the index ct m = -Again the constants, and i? m , are obtained by using the continuity 

conditions at the transition time rj eq between the radiation and matter era, /+ a d(?7eq) = Atm a t(?7eq) and /Lt' ad (?7eq) = 
A t mat(?7eq) , ) the full expressions are cumbersome and will not be presented. 

In Figure 4 we plot |/i mat |. We can see that a behavior like |/i ra d| is also obtained. As one may conclude, 
this happens because of the same reasons already discussed in the preceding section. For low frequencies if we 
have two perturbations |/ii(/3i)| and |A- 2 (^ 2 )|, with /3\ > /3 2 , then |/ii| > |/i 2 |. This condition inverts for the high 
frequency regime (for /?i > /3 2 then \hi\ < |1)- 
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FIG. 4: The amplitude of the GW in the matter era. a' 
conformal time for some selected values of /3 and a fixed 
function |/i ma t| for a high frequency, n = 10 4 . c) Modulus 
selected values of /3. 
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I Modulus of the mode function |/i mat | as a function of the 
comoving low frequency, n — 1. b) Modulus of the mode 
of the mode function |/i ma t| as a function of n and the same 
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FIG. 5: The power spectra of the GWs in the radiation and matter dominated eras. 


V. POWER AND ENERGY DENSITY SPECTRA 

Let us now discuss the power spectrum and the spectral energy density parameter (per logarithmic wave 
number interval) which can be written as (see Appendix for details): 


V{n,r)) 


32 G o 
- n 6 

7r 


IM’?)I 2 , 


(31) 


and 


fi gw (n,rj) = 7 ^^^(\K(v)\ 2 +n 2 \h n (r])\ 2 ). (32) 

In comparison to other alternative cosmologies, one advantage of our simple decaying vacuum models is that 
the above quantities can analytically be calculated for the radiation and matter dominated eras. Actually, this 
happens because the amplitude h(n 7 rj) = \ija and the corresponding solutions of f. i were explicitly obtained for 
each case (see previous section). 

In Figure 5 we show the plots of the power spectra, 'P ra d and P ma ti for radiation and matter dominated 
eras, respectively. From the first plot we see that TVad is almost flat until some transition frequency when 
begins to decrease. Similarly to what happens for the amplitudes, the power spectrum in this regime is slightly 
larger as the parameter /? increases. As should be expected, the decaying vacuum contributes to the creation 
of low-frequency gravitons with the corresponding spectrum remaining essentially flat. 

Nevertheless, after some transition frequency (corresponding ton ~ fewx 10 3 ), the waves are strongly damped 
and the associated power spectrum is no longer flat. This means that the decaying vacuum in this regime 
contributes more to increase the scale factor than to the production of gravitons. TAad decreases exponentially 
and the effect is even more pronounced for larger values of (3. Note also that at late times, the 'P-m&t for the 
matter-vacuum dominated phase presents the same high and low frequencies general properties of the radiation 
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FIG. 6: The energy density spectra of the GWs for the radiation and matter dominated eras. 


era. It starts with an almost flat spectrum and also decreases faster as long as the vacuum contribution is 
relatively larger (higher values of 0). 

In Figure 6 we display the energy density parameter as a function of the comoving wave number for the 
radiation and matter dominated eras. An interesting feature of the energy density spectrum in the radiation 
era, f2gw d \ is that if 0 > 0 then Hgw d \/3i) < &gw A \0) for all frequencies. Note also that figw^ grows as 
a power-law being weakly dependent on the value of 0 but a more strong dependence is obtained at the high 
frequency limit. 

The evolution of f2gw at ' follows a similar trends with some peculiarities at the high frequency limit. As in 
the radiation case if 0 > 0 then figw^^i) < Hgw at ^(/3 2 ) for all frequencies. In addition, for low frequencies 
the spectrum also grows as a power law linear being slightly lower for bigger values of /?. However, in the high 
frequency regime, the spectrum always decreases but varies differently as a function of the j3 parameter. In 
particular, for /3 ~ 0.2, the fall is very abrupt. Note also that for j3 = 0 it initially decreases and after remains 
almost flat for all modes. The basic reason for such a behavior is simple: when \h'\ 2 -C \h\ 2 the energy density 
spectrum is Sl gw oc n 5 \h\ 2r h~ 2 , where H = a'/a (see Eq. (B7)). This means that the contribution of the factor 

7~L~ 2 (which is not present in the power spectrum) largely determines the behavior of 


VI. FINAL COMMENTS 

In this paper, by assuming a spatially flat geometry, we have investigated the production of GWs for an 
interacting mixture of matter and vacuum in the context of GR. The dynamical A-term was described by a 
phenomenological law: A (H) = A& + 3/3H 2 , and a three stage description involving inflation, radiation and 
matter eras was adopted. It was also assumed that the vacuum decays only on the dominant component. 

For each cosmic era, we have determined the general expressions for the scale factor (in the conformal time), 
as well as the analytical solutions for the GW equation (see Figs. 1, 3 and 4). We notice that the mode 
function equations of the primordial GWs were derived and explicitly solved for each era. More interesting, 
the corresponding power spectra and the energy density parameter for the radiation and matter era were also 
obtained (see Figs. 5 and 6). Obviously, exact solutions are allowed in this framework due to the simplified 
phenomenological form adopted for the decaying A(I7)-term. 

In the present running vacuum model, the scale factor expands faster as long as the /3 parameter increases 
thereby affecting considerably the production and evolution of the GW modes for both regimes - low and high 
frequency limits (see Figs. 3 and 4). However, the most prominent feature is that the “potential” in the radiation 
era never vanishes (see Fig. 2), even when very small values of the /3 parameter are adopted. As a consequence, 
unlike models with no decaying vacuum (/3 = 0), the GWs in the radiation era can be adiabatically amplified 
(in the sense of Grishchuk [ 13]). This is the main result of the paper. It is different from the standard Parker’s 
result [46, 47] claiming that in relativistic cosmology there is no gravitationally induced quantum production of 
particles due the cosmic expansion in the the FRW radiation phase (for a similar result outside of GR see [48]. 
It is closely related with the fact that during the radiation phase a" ^ 0 in the presence of a decaying vacuum 
component. 

On the other hand, our expressions for the modulus of the GW amplitude \h\ shown us that the adiabatic 
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amplification (graviton production), is a low-frequency phenomenon even for this kind of A (H) cosmologies. 
This interesting and known feature is also maintained in the present context even considering that GW are 
produced in the radiation phase. Here as there, the basic problem is that at the high frequency regime the 
cosmic expansion dominates, and, therefore, the perturbations are dynamically suppressed in the course of the 
expansion. This behavior is also reproduced in the power spectrum due V ex n 3 \h\ 2 . However, in the case of the 
energy density parameter spectrum, H gw oc 'H~ 2 n 2 V , the behavior is different due to the contribution of TL~ 2 . 

Finally, we stress that the present work was based on a very simple decaying vacuum cosmology. It played 
the role of a toy model for obtaining analytically some basic information, like the production of GWs during 
the radiation phase. Its consequences on the B-modes polarization of CMB anisotropies it will discussed in a 
forthcoming communication. As it appears, the model can be thought as a starting point for the investigation 
of more complex and rich decaying vacuum cosmologies, like the one proposed in Refs. [25, 49], and recently, 
investigated in a more enlarged framework in Refs. [29, 30]. Although less analytical regarding the calculations 
involving GWs, such models deserve a closer scrutiny since they furnish a complete cosmological history. 


Appendix A: Quantized tensor perturbations 

The generating mechanism of the primordial GWs is believed to have a quantum mechanical origin. The basic 
idea is that quantum fluctuations of the vacuum state in the early universe were stretched to macroscopic scales 
due the cosmic inflation thereby originating the present observable primordial GW spectrum. The standard 
quantization procedure is based on a semi-classical approach where the perturbations are quantized but the 
gravitational background evolves classically (for details see the review of Giovannini [50]). 

When the perturbations are quantized on a classical background, the functions c n e c n in (18) are promoted 
to be quantum creation and annihilation operators which satisfy equal time commutation relations (H = 1 ) 

1^,^] = W 3 (n-m), c„|0) n = 0. (Al) 

The vacuum state |0) n is defined for a given time 77 and mode n. In general does not exist an unique 
vacuum state because if for some time 7ji the vacuum state is c n (r]i) | 0 j) n = 0 , for another time rjf it changes to 
c n (b/)|0/)n 7 ^ 0. For an expanding background, the vacuum state at time r]i is different from the vacuum state 
at rjf. This so-called vacuum state ambiguity has some interesting consequences, the most prominent of them 
being the gravitationally-induced particle production (‘Grishchuk particles’). To be more precise, by assuming 
that there is no particles at time rji, that is, W| 0 i) = 0 , where W is the number operator acting in the vacuum 
state [Nt = c\ci (the polarization and wave number index were suppressed for simplicity)]. Later on, at a time 
T)f we have Nf\0f) = 0, but the vacuum state ambiguity gives iV/|0j) ^ 0. Consequently the vacuum state |0j) 
contains “/” particles and vice versa. This is the physical foundation for the creation of ‘Grishchuk gravitons’ 
induced by the cosmic expansion. 


Appendix B: Power and energy density spectra 

With the perturbations already quantized, the important physical observables are readily calculated. The 
(dimensionless) power spectrum, that is, the quadratic mean value of the amplitude of the perturbations, can 
be defined as: 


V(n,ri) 


d(0|hi J (77,x)/W(77,x)|0) 

dlnn 


Now, by inserting (Al) into (18) one finds 


(Bl) 


(0|/i.„-(77,x)/i lJ (7y,x)|0) 


— [°° n 3 \h n (ri)\ 2 dhin. 


7r 


o 


(B2) 
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Note that the equality holds because the modulus of the modes for different polarizations are equal, \h\ 2 = 

X 

\h \ 2 = \h\ 2 . Finally we have: 


V{n,r]) 


32 G 

7 T 


n 3 \h n (v)\ 2 . 


Another important quantity is the energy spectrum defined in the following manner: 


(B3) 


v) 


1 d(0\p sw (r])\0) 
Pcrit dlnn 


(B4) 


which represents the GW energy density (p gw ) per logarithmic wave number interval, in units of the critical 
density p C nt(il) = 3H 2 (^)/8nG. The GW density is 


^ i K ) 2 + ( vft ,,) 2 

Psw _ 0 _ 64t rG a 2 

whose vacuum expectation value in the vacuum state reads: 

<°l P ‘ M = J 0 2^-- H- 

while the energy spectrum (B4) becomes 


f2g W (ri) tf) 


8nG n 3 

mHrf) 2^2 


(|/ l ;(, ? )| 2 +n 2 |/ l „(7 ? )| 2 ). 


Both quantities are the most important primordial GW observables. 


(B5) 


(B6) 


(B7) 
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